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Abstract 

The hitherto unobserved purely leptonic decays — ► t^v t and B^ — ► \i Vn are of much 
interest at current and future runs of the e + e~ B factories. Such decays are sensitive 
to charged Higgs bosons (H^) at the tree-level and provide important constraints on 
tan/3/m#±. We include the large corrections to the H^ub coupling induced by virtual 
SUSY effects and show that the bounds on tan/3/m^± can be significantly weakened or 
strengthened. 
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Figure 1: The decay — > i^vi is mediated by annihilation into a virtual W (and in extensions 
of the SM also by a charged Higgs). 

1 Introduction 

The purely leptonic decays of the charged B mesons, B — > r^v-y and B ± — > y u^, are of much 
interest, both in the Standard Model (SM) and in the context of models beyond the SM. The 
vast majority of B ± mesons (i.e. a bound state of bu or bu) decay via the subprocess b — > qX, 
where q = c,u,s,d (the latter two via penguin operators) and X are either quarks or leptons. 

A purely leptonic final state can occur via the annihilation process B^ 1 — > W* — > i ! i/j, 
although with a very suppressed rate (fig. Q). Annihilation to y is the dominant decay 
mechanism for the lighter charged mesons 7r ± and K , while for Df the purely leptonic decays 
t^Vt and y^v^ have been observed with subdominant branching ratios (BRs) in the range 
10 -2 — > 10 -3 pp. No such purely leptonic decay has been observed for B ± mesons. Although 
their search is problematic at the hadronic B factories (e.g. Tevatron Run II) due to the missing 
energy of the z/, upper limits on their BRs were obtained at e + e~ colliders by CLEO 0111! and 
the LEP collaborations |4 j. BELLE and BABAR [H] are already producing more stringent 
limits, and expect to be sensitive to the SM rate by 2005/2006. Observation of such decays 
would offer a direct measurement of the B ± decay constant 

These decays are also sensitive to many well motivated extensions of the SM jjj. Charged 
Higgs bosons (H^) appear in any model with two Higgs SU{2) x U{1) doublets (which includes 
all supersymmetric models) and/or with Higgs triplets (which can provide neutrino mass) 
In contrast to many other rare decays (e.g. b —>■ s'j ) which are influenced by new physics at 
the one loop level, the decays i? ± — > l^vi are sensitive at tree-level to H ± and can be enhanced 
up to the current experimental limits [7j. From these upper limits very useful constraints 
on the parameter tan/?/m^± can be obtained (where tan/5 = t^/fi, and i>j are the vacuum 
expectation values of the Higgs doublets), which are comparable (and often superior) to those 
from processes like B — > Xr^v [5] and t — > H^b. The SM predictions for B ± — > l^vi suffer 
from a sizeable uncertainty in Jb, whose value can only be calculated using non-perturbative 
techniques such as lattice QCD or inferred from other processes (e.g. B°B° mixing). The 
current errors in such calculations induce an uncertainty in the SM BRs of around 30% [TUj . It 
is likely that this sizeable error together with the small rates has dissuaded explicit calculations 
of higher order corrections, which would naively be expected to be at the 10% level or less. 
However, we turn out attention to B ± —>■ r^v T and B^ —>■ y for three reasons. 

(i) Large corrections to the H^ub coupling are possible in supersymmetric (SUSY) models for 
large tan /3 [T2"lll3lll4| . Such corrections should significantly affect the BRs for B ± — > t ± v t 
and B ± — > y^v^ and the derived constraints on tan (3/m H ±. 
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Table 1: SM predictions and current experimental limits from various machines. 



(ii) Due to the successful running of the B factories, the SM rate for both decays will be 
within experimental sensitivity by 2005/2006. Higher Luminosity upgrades are now being 
discussed [TS] which would offer an ideal environment for precision studies of these decays. 

(iii) Improvements in lattice techniques resulting in more precise predictions for fs jTE] are 
expected in the next few years. Data from the recently approved CLEO-c jTTj will play 
an important role in refining the techniques. 

The advances in (ii) and (iii) will facilitate the search for new physics in 5 ± — > r ± i/ T and 
B^ — > [i Up, and leads us to consider the higher order corrections to these decays in (i). 

Our work is organised as follows. In Section 2 we briefly review the tree-level contribution 
of H ± and summarise the current experimental status. In Section 3 we introduce the SUSY 
corrections while Section 4 presents our numerical results. Finally, Section 5 contains our 
conclusions. 



2 B ± — ► r ± i/ T and B ± — ► 



The and H ± induce effective four-fermion interactions [7] for bu —> I v% and bu — > l + vi. 
Below we write the Lagrangian for bu — > l~T>i induced by W~ and H~ exchange: 

([«7^(l - 7fe)6] [W-TsH -Ri[u(l + K)b] KI-75H) (1) 

with 

Ri = tan 2 (5 (iribmi / m 2 H ±) (2) 

Here we only keep the dominant mb tan /3 part of the H ± ub coupling. The two quarks need 
to be at zero separation in order to annihilate and this probability is contained in the decay 
constant which is defined as: 

< 0\uj^ 5 b\B- >= if BV %- < 0|U 7 5 6|5- >= -if B (m%/m b ) (3) 

The decays proceed via the axial-vector part of the coupling and via the pseudoscalar part 
of the H ± coupling. Other couplings give a vanishing contribution. The decay amplitude is 
therefore: 

M = -^JfubfB [mi - R^ml/rrib)] 7(1 - (4) 
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The tree-level partial width is given by [7]: 
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(5) 



where 



r H = [1 - tan 2 f3(m 2 B /m 2 H ±)] 2 



(6) 



The overall factor of mf arises from the helicity suppression of the W ± contribution, and the 
Yukawa coupling of H I Vi. 1 The dependence on m& cancels out in the tree-level approx- 
imation. This is because mj in Ri originates from the b quark Yukawa coupling (y&) which 
at lowest order is simply related to m;, by yj, = V^rrib/vi. Therefore m& in the numerator 
in M. cancels with m;, in the denominator. At the 1-loop level the simple relationship be- 
tween yi, and the running quark mass m& is modified, resulting in large corrections which will 
be discussed in Section 3. One can see from the expression for r# that the H ± contribution 
interferes destructively with that of W^. Enhancement of these decays (r# > 1) requires 
tan(5 /m B ± > 0.27 GeV -1 , while exact cancellation (r# = 0) of the W ± and H ± contributions 
occurs for tan(3/m H ± « 0.19 GeV -1 (= 1/M B ). 

Table displays the SM predictions for BR(i? ± — > l^vi) along with the current experimental 
limits from various machines. We take f B = 200 MeV and \V u b\ = 0.0035. The BRs are 
hierarchical due to the aforementioned mf suppression. BR(5 ± — > r ± z/ r ) is much larger than 
BR(£> ± — > /i ± z/^), but both decays are within range of data samples with C(10 8 ) mesons. 
BR(5 ± — > e u e ) is too small for any current or planned experiment. The detection efficiency 
for B ± — > [J, Up is far superior to that for B^ — > t^Vt, which compensates for the smaller BR 
of the muonic channel. Hence these two decays are competitive with each other for probing 
new physics and/or for offering measurements of the decay constant f B . The upper limits in 
Table 1 on 5 ± — > r v T (from BABAR) and B^ — > \l (from BELLE) constrain r B < 4.7 
and 17.4 respectively. For f B = 200 MeV this gives a bound on tan. /3 /run < 0.34 ±0.02 GeV -1 
from B^ —>■ r^Ur, with weaker limits from B ± — > y^Py,. Data samples of order 4 x 10 8 B ± 
mesons, which are expected by 2006, should allow sensitivity to the SM rate (r# = 1). BABAR 
simulations give yields of 17 and 8 events (after all cuts) for the t^Uj. and \i channels 
respectively (13], which would be enough to make a first measurement of the BRs, or, in the 
case of non-observation, would restrict tan/?/mjj± < 0.27 GeV -1 . 

3 Higher order corrections 

Although the tree-level expression contains a sizeable error from the uncertainty in f B V u b we 
feel that including higher order corrections to the rates for B^ — > r^Ur and B ± — > \i 
is well motivated, since current experiments will reach the sensitivity of the SM prediction 
in the next few years. In addition higher luminosity upgrades of both PEP-II and KEKB 
are being discussed, which would provide data samples of 10 9 — > 10 10 B ± mesons and allow 
precision measurements of these decays |15j . Importantly, the contribution can already 
saturate the current experimental limits (as shown in Section Ej), and so one or both of t^v t 

1 Non-TO/-suppressed purely leptonic decays are possible in other models, e.g. i?-parity violating scenarios 
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and y^v^ might be observed in the near future. The constraints on tan/5/m^± derived from 
these decays are among the most model independent constraints on this important ratio. We 
stress that other rare decays (e.g. B — > X s j, B^ s — > H + H~) can also give strong constraints on 
tan/3/m H ± [13111112011211 • 

The observed inclusive decay B — > X c r v [22] is also sensitive to if at tree- level and 
provides competitive, mostly model independent limits of tan /3/m H ± < 0.5. The exclusive 
decay i? ± — > D°r ± z/ will be pursued at the B factories |23j . Constraints on tan f}/mH± from 
the unobserved top quark decay t — > H^b at the Tevatron Run I are inferior to those from 
B ± -> r±f and 5 ± -> X c r ± z/, but will improve significantly with 2ft) 1 expected during Run 
II [21]. Given the importance of these constraints on tan/3/m#±, we feel that higher order 
corrections to the rates for BR(S ± — > t^v) and BR(S ± — >• [V s1 v) should be considered. We 
believe that a closer analysis of these decays is timely due to the anticipated improvements in 
lattice calculations of Jb in next few years [To] . 

In this study we will only consider the sizeable corrections to the coupling H ± uidj which 
occur at large tan/3 in SUSY models. The decays of interest to us depend on the coupling 
H^ub (I = 1, J = 3) and the corrections are comprised of the following [2*5] : 

(i) Enhanced corrections to H ± ub vertex diagrams 

(ii) Large corrections to the b quark Yukawa coupling yt, [T21["""""] (known as the HRS effect) 

(iii) Large corrections to V u b [T"I], where V u b in the MSSM Lagrangian differs from the exper- 
imentally measured 

The above corrections greatly affect other processes involving the H^uidj coupling such 
as b — > S7 [26J and t — ► H^b [2Zj (see [2H] for a review). The effects of these corrections on 
B ± — > XcT^u were covered in [21] but the decays B ± — > f^v T and B ± — > of interest to 

us were overlooked. It has been shown [22] that the combined effect of these corrections can be 
neatly encoded in an effective Lagrangian for the H^ujdj coupling by multiplying the usual 
Lagrangian by a factor 1/(1 + ejtan/3). For our case of the decay of a B u , J = 1 and the 
corresponding ij is called e in |25j. 

This correction appears in Ri (eq"2J) and modifies the scaling factor of the SM rate to: 

H \ 1 + e tan (3 m 2 H± ) ^ 

Analogous corrections of type (i) and (ii) are also present for the H ± r ± v T and H^/jtv^ 
couplings, but are known to be considerably smaller since there is no SUSY QCD contribution 
(only SUSY electroweak contributions). Thus we will neglect them in the present study. These 
corrections will be considered in a future work and are of more interest for high luminosity 
runs of the B factories with data samples of 10 10 B^ mesons. Note that both fg and the large 
corrections to the H ± ub coupling cancel out in the ratio i? MT defined by: 



BR{B ± -> r 



R » T ~ BR(B ± - p±ul) (8) 

This quantity is predicted to be 0.8 yl/y^, which in the SM is essentially equivalent to 0.8 m^/m^. 
However, R^ T would be affected by the HRS corrections to y T and the vertex corrections to 
z vi. 
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Figure 2: Excluded regions in the tan (3 — rriH± plane for different values of eo and fs from the 
experimental limit BR(5 + — > tu) < 4.1 • 10~ 4 . The corresponding plot for the n case looks 
very similar but gives slightly weaker constraints. 



4 Numerical Analysis 

In this section we discuss the numerical impact of the modification of the H^uh vertex on 
the combined bounds on tan/? and m#±. We take the modification parameter e as an input 
parameter, varying it in the interval [—0.01,0.01]. This is the range found in [251 by explicit 
calculation and a scan over reasonable values of the Minimum Flavour Violation Model param- 
eters. 

In order to show the region in the tan/3 — m^± plane excluded by the experimental upper 
limits on BR(5 ± — > t^v t \ in FigJ2]we plot tan/3 (y-axis) against (x-axis). For a given 

set of input parameters, the region above the corresponding curve is excluded. 

From BR(i? ± — > r u T ) < 4.1 ■ 10~ 4 , without the radiative corrections, (eo — 0) the excluded 
region is simply bounded by a straight line of gradient R = 0.34 (0.36, 0.32) [GeV -1 ] for 
fs = 200 (170, 230) MeV. With the inclusion of radiative corrections at the level of e"o = —0.01 
(0.01), the bounds change to the lower (upper) set of curved lines. It can be seen that even 
with the present rather large uncertainty of Jb of ±15%, the effect of the radiative corrections 
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is much larger than that of the uncertainty of the decay constant. 

With the central value of fs = 200 MeV, m H ± = 120 GeV would allow tan (3 of up to 40 
when the radiative corrections are neglected, strengthening to tan/5 < 33 for e = —0.01 and 
weakening to tan/3 < 49 for e = 0.01. The corresponding bounds from BR(5 ± — > fi u^) < 
6.5 • 10~ 6 are slightly weaker: tan/? < 51 for e = and 40 (66) for e = -0.01 (0.01) (all for 
fs = 200 MeV, m H ± = 120 GeV). Variations of fs have a much smaller effect. For the r case 
and with m H ± = 120 GeV, the bound on tan (3 changes from 40.5 (200 MeV) to 42.8 (38.6) for 
f B = 170 (230) MeV. 

5 Conclusions 

The unobserved purely leptonic decays B ± — > t ± u t and B ± — > /i ± ^ provide competitive 
constraints on the ratio tan/?/m#±. We have studied the effect of SUSY corrections at large 
tan/? and showed that the exact bounds on tan/3/m#± can be strengthened or weakened 
depending on the value of e , the correction factor for the H ± ub vertex due to SUSY loop 
corrections. The effect on the allowed region in the (tan/3, m#±) plane exceeds ±20% for 
eo = ±0.01, which is considerably larger than the analogous uncertainty of ±5% caused by 
varying the decay constant fs by ±30 MeV. 
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